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Abstract 

Two approaches were taken to make convenient spread sheet calculations of elastic constants 
from resonance data and the tables in ASTM C1259 and E1876: (1) polynomials were fit to the 
tables; (2) an automated spread sheet interpolation routine was generated. To compare the 
approaches, the resonant frequencies of circular plates made of glass, hardened maraging steel, 
alpha silicon carbide, silicon nitride, tungsten carbide, tape cast NiO-YSZ, and zinc selenide 
were measured. The elastic constants, as calculated via the polynomials and linear interpolation 
of the tabular data in ASTM Cl 259 and El 876, were found comparable for engineering 
purposes, with the differences typically being less than 0.5 percent. Calculation of additional v 
values at t/R between 0 and 0.2 would allow better curve fits. This is not necessary for common 
engineering purposes, however, it might benefit the testing of emerging thin structures such as 
fuel cell electrolytes, gas conversion membranes, and coatings when Poisson’s ratio is less than 
0.15 and high precision is needed. 


Introduction 

The elastic constants of isotopic materials can be measured by using a variety of static and 
dynamic techniques such as strain gages, impulse excitation, or resonant ultrasound 
spectroscopy, and a recent comparison of the techniques has been made (ref. 1). Dynamic 
techniques have the advantage of simple test specimen geometry, good precision, and 
applicability over a wide range of temperatures. One convenient method is the impulse excitation 
technique in which a test specimen such as a rectangular beam or circular plate is excited by a 
brief mechanical impulse with a tool (see fig. 1) and allowed to vibrate naturally. The specimen 
geometry and support and impact locations are chosen to induce flexural or torsional modes of 
vibration. The vibrations are detected with a non-contacting microphone or contact 
accelerometer, and are analyzed by a signal processor to determine the fundamental resonant 
frequency. The fundamental frequency is used with the specimen geometry and mass to estimate 
dynamic elastic properties such as Young’s modulus, Poisson’s ratio, and shear modulus. 

This method is nondestructive, sensitive, requires relatively little material, and can be used to 
test specimens designed for other purposes such a uniaxial or biaxial flexural strength 
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measurement (refs. 2 to 4). It also is useful for monitoring changes in newly developed materials 
that contain large amounts of porosity and undergo chemical changes such as reduction or 
oxidation during testing and use. Another advantage of the technique is the availability of 
American Society for Testing and Materials (ASTM) standards that give detailed guidance on 
the technique (refs. 5 and 6). 

For a circular plate, the resultant first and second natural vibrations corresponding to torsion 
and flexure are measured and used to independently estimate Young’s Modulus. The values are 
then used to calculate an average Young’s modulus with the relation (ref. 5) 

E = {E, + E 2 )/2 (1) 


where: 


E,= 


37.6991 f 2 P 2 m{l-v 2 ) 


K]t 3 


and E 2 = 


37.6991 f;D 2 m(l-v 2 ) 


Kit 3 


( 2 ) 


with fi and /2 being the first and second natural frequencies, D the disc diameter, t the thickness, 
m the mass, u Poisson’s ratio, and K } and K 2 the first and second natural geometric factors. The 
shear modulus, G, is calculated from the common relation 


G = 


E 

2(1 + v)' 


( 3 ) 


Currently, the standards ASTM C1259 and El 876 (refs. 5 and 6) use linear interpolation of 
data compiled in tables to determine the constant v as a function of t/R and fi/f 2 where R is the 
disc radius. The constants Kj and K 2 are then determined with linear interpolation as a function of 
t/R and v. This is inconvenient as three manual interpolations are required per constant, and the 
accuracy is reduced because the functions are nonlinear and the interpolations of Ki and K 2 
depend on those of v. Although extremely accurate values of E and v are not needed for general 
engineering purposes, accurate constants are particularly useful in verification of test rigs and 
models. Further, rapid, convenient calculations are beneficial if elastic constants are used to 
nondestructively monitor material changes. 

In order to allow accurate, convenient spread sheet calculations of elastic constants from 
plate data, two approaches were taken: (1) polynomials were fit to the data tables A 1.1 to A 1.5 of 
ASTM C1259 and El 876, which are identical and based on the data of Glandus and Martincek 
(refs. 7 and 8); and (2) an automated bilinear spread sheet interpolation routine was implemented 
(see the app.). The resultant polynomials exhibited residuals of less than 0.5 percent for A; and 
K 2 , and less than 1.5 percent for v, and are sufficient for most engineering purposes. 

Comparisons of calculations made with linear interpolation and the polynomials were made for 
glass, steel, silicon carbide, silicon nitride, tungsten carbide, tape cast NiO-YSZ and zinc 
selenide, and indicate similar results. 
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Polynomial Functions 


During fitting of the constants, an attempt was made to minimize the residuals while using a 
relatively simple polynomial function. This was possible for both K t and K 2 , however reasonable 
fits could not easily be made for v. The most difficulty in fitting v was encountered for values of 
f 2 /fi = 1.350, for which v varies by a factor of ~6 and changes from curvilinear behavior to very 
linear behavior at small v as shown in figure 2. Thus three separate functions were generated to 
describe the range of v given in table A 1.1 of C1259. In addition, a wide range fit to the data of 
Glandus (ref. 7) that did not include values of v for f 2 /f 2 = 1.350 was generated. Other regions of 
the table that presented difficulty were additional columns containing several constant values of 
v, as with th ef 2 /fi = 1.350 data, and columns containing linearly increasing values of v as a 
function of t/R. Such regions of constant v or linearly changing v did not follow a smooth trend 
as might be expected for a physical phenomenon. This may have resulted from the numerical 
accuracy available when Martincek did the calculations and the use of a nomogram to present the 
data, and the fact that many of the columns and rows in table A 1.1 were generated by 
interpolation of the data of Glandus (ref. 7), rather than from direct calculation. Refinement of 
the calculations for v as a function of t/R and fVfi might allow better curve fits over a wider 
range. 

An example of constant and linearly varying v is shown in figure 2. Unfortunately, these 
complications tend to occur at small t/R ratios, which are commonly used in the mechanical 
testing of structural ceramics such as membranes, fuel cell elements, and coatings that are being 
manufactured in very thin sections (t/R < 0.05) and are required to carry thermal and structural 
loads. Fortunately, the complications occur for v values (< 0.15) not usually exhibited by 
ceramics (typical v = 0. 16 to 0.30) and the curves are relatively flat. For measurements of thin 
structures exhibiting low Poisson’s ratio, improved values of v between t/R = 0 and 0.2 might be 
beneficial. 


For 0 <t/R< 0.500, v can be described for 1.375 <f 2 /fi < 1.525, 1.525 <f 2 /fi <1.725, and 
1.725 <f 2 /fi < 1.900 by the expression: 


v = a + b\ 


f t ' 


\Rj 

+ h 


+ c 


( I?" 

\fij 


+ d\ 


f t^ 2 


\Rj 


+ e\ 


\fi J 
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ro 2 
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r A'" 

\fi j 


+g 




\Rj 


( 4 ) 


where the constants a through j are given in table 1 . The fits are within 1 percent of the data in 
table Al.l C1259. For the wider range of 0.00 < t/R < 0.500, and 1.400 < f 2 /fi < 1.900, v can be 
estimated within 1 .3 percent of the data of Glandus (ref. 7) by using the constants in the last row 
of table 1 with equation (4). The v values corresponding to f 2 /fi = 1.350 were not included in the 
fit because they decreased the quality of fit and are not commonly encountered in experimental 
measurements of ceramics. Examples of fits to the data are shown in figures 3 and 4. 

For fitting of K/ and K 2 data, minor difficulties occurred around some data points. Although 
Glandus added many intermediate data points to the range of Kj and K 2 given by Martincek, he 
appears to have re-used Martincek’s data directly without adding precision: the data in Glandus’s 
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tables are given to three places; however, many of the data points are identical to those of 
Martincek which were originally reported to two places. Evidently, Glandus simply added a third 
decimal (i.e., a zero) to Martincek’s data rather than recalculating those points to a higher 
precision. 


TABLE 1.— CONSTANTS FOR ESTIMATION OF V FROM EQUATION (4) 



fjfi range 

Coefficient 

1.375-1.525 

1.525-1.725 

1.725-1.900 

1.400-1.900 

a 

0.7161 

-4.3343 

-8.0933 

-4.30526 

b 

-0.6502 

2.1017 

3.4452 

0.38999 

c 

-4.1377 

5.6718 

12.2471 

6.10124 

d 

1.0459 

0.5107 

-2.2461 

0.19582 

e 

4.2095 

-2.2025 

-6.1198 

-2.76353 

f 

0.6056 

-2.5965 

-3.2070 

-0.47553 

% 

-0.5148 

-0.5109 

-0.0872 

-0.62008 

h 

-1.1313 

0.2792 

1.0725 

0.45518 

i 

-0.1247 

0.7924 

0.7117 

0.11816 

j 

-0.3197 

0.0283 

1.4525 

0.31933 


For 0 <t/R< 0.500 and 0 < v< 0.500, Kj and K 2 can be described by the expression: 


( * h 


a + b\ 


K. = 


\R; 


+ c 


ft } 2 

kRj 


+ dv + ev 2 + fv 3 


1 + g 


f -1 


+ h 


f -1 


+ iv + jv 


( 5 ) 


where the constants a through j are given in table 2. The fits are within 0.5 percent of the data in 
table A1.2 and A1.4 of C1259. For the narrow range of 0.100 < t/R < 0.200 and 0.14 < v< 0.34 
given in tables Al. 3 and A1.5, Kj and K 2 can be described by the expression: 


K , = a + b\ 


\Rj 


+ cv + d 


\Rj 


+ ev 2 + / 


\Rj 


V + g\ 


OV 

\Ry 



(t^ 

2 

ro 

+ i 


V- +J 



kRj 


UJ 


( 6 ) 


The fits are within 0. 1 percent of the data. Examples of fits to the data are shown in 
figures 5 to 8. 
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TABLE 2.— CONSTANTS FOR EQUATIONS (5) AND (6) 



0 <t/R< 0.500 and 0 < v< 0.500 

0.100 <t/R< 0.200 and 0.14 < v<0.34 

Coefficient 

Kj 

Ki 

K, 

k 2 

a 

6.1528 

8.2584 

6.2305 

8.3901 

b 

0.5085 

-2.3794 

-1.7997 

-2.9434 

c 

0.8825 

6.6164 

-2.4825 

2.0928 

d 

-3.9037 

2.3540 

3.6354 

5.0895 

e 

-5.4739 

-9.1109 

-0.3022 

1.9810 

f 

3.9297 

-0.9923 

3.6488 

2.7693 

s 

0.1364 

-0.2142 

-9.3417 

-18.3125 

h 

0.9120 

1.7467 

-2.2626 

-3.0307 

i 

-0.3628 

-0.0475 

-3.3588 

-3.1336 

j 

-0.3612 

-0.9262 

-10.7968 

-15.1250 


Experimental Results 

In order to directly compare test results using linear interpolation and the polynomials, 
circular plates machined in accordance with the procedures of ASTM Cl 161 and C1499 
(refs. 1 and 2) were tested. The dimensions and mass of all the test specimens were measured 
with a resolution of 0.001 mm and 0.0001 g, respectively. The plates were typically 50. 1 mm in 
diameter and ranged in thickness from 1.54 to 4.53 mm ( t/R = 0.06 to 0.18), and were intended 
for strength testing. The materials included glass, hardened maraging steel, alpha silicon carbide, 
silicon nitride, tungsten carbide, NiO-Yttria-stabilized zirconia (YSZ), and zinc selenide. The 
75mol%NiO-YSZ was an unreduced anode material with ~32 percent porosity. It was 
manufactured by tape casting and sintered at 1200 °C for 2 hours as part of a fuel cell program 
being conducted at NASA Glenn Research Center (ref. 9). The zinc selenide was optical grade 
material made by chemical vapor deposition and intended for use on the International Space 
Station (ref. 10). 

The test specimens were supported on the torsional or flexural nodal lines, as necessary, with 
foam rubber supports as shown in figure 1 and lightly impacted with a steel tipped hammer. The 
results are summarized in tables 3 and 4, and indicate that the differences in Poisson’s ratio and 
Young’s modulus as estimated from interpolation and the polynomials are usually less than 0.5 
percent. 


Conclusion 

Polynomial fits and linear interpolation of the data tables of ASTM C1259 and El 876 
produce comparable estimates of Poisson’s ratio and Young’s modulus for engineering purposes 
For glass, hardened maraging steel, alpha silicon carbide, silicon nitride, tungsten carbide, NiO- 
YSZ, and zinc selenide the differences as calculated from polynomials and interpolation were 
typically less than 0.5 percent. Calculation of additional v values at t/R between 0 and 0.2 would 
allow better curve fits and might benefit the testing of emerging thin structures such as fuel cell 
electrolytes, gas conversion membranes, and coatings when Poisson’s ratio is less than 0.15 and 
high precision is needed. However, the current values are sufficient for common engineering 
purposes. 


NAS A/TM— 2005-2 1 3843 


5 





Appendix 


The following Microsoft® Excel macro 1 was used to linearly interpolate values ofu, Kj, and 

K 2 . The interpolated values were used with the physical data in tables 3 and 4 to estimate E t and 

E 2 which were than averaged to find the Dynamic Young’s Modulus (E). 

’USER DEFINED FUNCTION 

’XY interpolate 

’Interpolates values from a named 2D table for row(x) and column(y) variables. 

’The named table should include the row and column variables. 

'Cell (1,1) should not contain a numeric value or formula. 

'Row variables should increase to the right. Column variables should increase downwards. 

Static Function XYinterpolate(xyarray As Variant, x, y As Single) As Single 
Dim nl, ml as integer 

Dim xl, x2, yl, y2, Rylxl, Rylx2, Rylxlx2, Ry2xl, Ry2x2, Ry2xlx2 As Single 

xl = Application.HLookup(x, xyarray, 1) 

nl = Application.Match(xl, xyarray.Rows(l), 0) 

x2 = xyarray.Cells(l, nl+1). Value 

yl = Application. VLookup(y, xyarray, 1) 

ml = Application.Match(yl, xyarray.Columns(l), 0) 

y2 = xyarray. Cells(ml+1, l).Value 

Rylxl = xyarray. Cells(ml, nl) 

Rylx2 = xyarray. Cells(ml, nl+1) 

Rylxlx2 = (x-xl)/(x2-xl)*(Rylx2-Rylxl)+Rylxl 
Ry2xl = xyarray.Cells(ml+l, nl) 

Ry2x2 = xyarray. Cells(m 1+1, nl+1) 

Ry2xlx2 = (x-x 1 )/ (x2-x 1 ) * (Ry2x2-Ry2x 1 )+Ry2x 1 
XYinterpolate = (y-yl)/(y2-yl)*(Ry2xlx2-Rylxlx2)+Rylxlx2 

End Function 

^his macro was courtesy of Dennis Kirk Engineering, denniskb@ozemail.com.au. 
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TABLE 3.— COMPARISON OF E AND V VALUES FROM INTERPOLATION AND NARROW RANGE CURVE FITS. VALUES AFTER ± ARE ONE 

STANDARD DEVIATION. PERCENT DIFFERENCE IS FROM THE INTERPOLATED VALUE 
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TABLE 4.-COMP ARISON OF E AND V VALUES FROM INTERPOLATION AND WIDE RANGE CURVE FIT. VALUES AFTER ± ARE ONE 
STANDARD DEVIATION. PERCENT DIFFERENCE IS FROM THE INTERPOLATED VALUE 
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Figure 1 . — Schematic of support setup for Impulse Excitation Technique: (a) torsional mode 
of vibration and (b) flexural mode of vibration. 
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Figure 2. — V as a function of t/R for regions exhibiting linear variation or no 
variation of V with t/R. The data is from table A1.1 of ASTM Cl 259 (ref. 4). 



0.00 0.10 0.20 0.30 0.40 0.50 

t/R 


Figure 3. — V as a function t/R for various f^f i- Solid lines are equation 4, solid symbols 
are data points from table A1 .1 of ASTM Cl 259 (ref. 4) and references 6 and 7. 
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t/R 

Figure 4. — V as a function t/R for various f^f i . Solid lines are equation (4) wide 
range fit, solid symbols are data points from table A1.1 of ASTM Cl 259 (ref. 4) 
and references 6 and 7. 



Figure 5. — K - 1 as a function t/R for various V. Solid lines are equation (5), solid symbols 
are data points from table A1 .2 of ASTM Cl 259 (ref. 4) and references 6 and 7. 
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Figure 6. — K- \ as a function tIR for various V. Solid lines are equation (6), solid symbols 
are data points from table A1.3 of ASTM Cl 259 (ref. 4) and references 6 and 7. 



tIR 

Figure 7. — K 2 as a function t/R for various V. Solid lines are equation (5), solid symbols 
are data points from table A1 .4 of ASTM Cl 259 (ref. 4) and references 6 and 7. 
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Figure 8. — K 2 as a function t/R for various V. Solid lines are equation (6), solid symbols 
are data points from table A1.5 of ASTM Cl 259 (ref. 4) and references 6 and 7. 
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